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Abstract 

We consider decay processes of scalar-field condensation in the framework of well- 
established quantum field theory. We postulate that the quantum state corresponding 
to the scalar-field condensation is so-called coherent state with discussing the validity 
of such a treatment. We show that, by using the unitarity relation of the scattering 
matrix, decay rate of the coherent state is systematically calculated. We apply our 
procedure to derive explicit formulae of decay rates for two cases: (i) we study the 
case where the scalar condensation decays into a pair of scalar particles and show that 
our formalism reproduces the results obtained from the parametric-resonance analysis, 
and (ii) we also calculate the decay rate when the coherent state decays via anomaly. 



1 Introduction 



Scalar-field condensations play very crucial roles in various places in cosmology. Probably 
the most important example is the infiaton field which is necessary for infiation [HE]- In 
particular, in the slow-roll inflation models [3l H], energy density of the inflaton conden- 
sation provides the energy density to realize (quasi-) de-Sitter universe during inflation. 
After inflation, the inflaton oscillates around the minimum of the potential and decays into 
standard-model particles to reheat the universe realizing hot big-bang cosmology. This class 
of inflation model not only solves flatness and horizon problems in cosmology but also pro- 
vides density fluctuation consistent with the Wilkinson Microwave Anisotropy Probe data [5] . 
In addition, it has been pointed out that the density fluctuation may arise from late-decaying 
scalar fleld other than the inflaton, which is called "curvaton" [HI [71 E]. Other important 
example is the Affleck-Dine fleld for baryogenesis [9J . In low-energy supersymmetric models, 
there exist scalar flelds, i.e., scalar partners of quarks, which have baryon number. If some 
of those flelds acquire non-vanishing amplitudes in the early universe, non- vanishing baryon 
number may be imprinted into the motion of the scalar-quark condensations due to baryon- 
number violating operators at an ultra-high energy scale. Such a scenario is one of the most 
attractive scenario to generate large enough baryon asymmetry of the universe. 

All of these exotic scalar-fleld condensations (i.e, inflaton, curvaton, Affleck-Dine fleld, 
and so on) oscillate around the minimum of the potential at some stage of the evolution 
of the universe, and eventually decay into standard-model particles for the cosmological 
history consistent with observations. Thus, it is important to understand how the scalar- 
fleld condensation decays from the view point of the quantum fleld theory. 

The main concern of this paper is to discuss how the oscillating scalar-fleld condensation 
decays into other states. Around the minimum, the potential of the scalar fleld ip is well 
estimated by a parabolic one 

V = \m'y. (1.1) 

Neglecting the effect of the cosmic expansion for simplicity, the solution to the classical 
equation of motion is given by 

if = A^coiim^t. (1.2) 

In this case, the energy density of the condensation is given by = ^m'^A'^. One should 
understand how the energy density stored in the oscillation of ip is converted to the that 
of radiation. In the simplest approach, the decay rate of the scalar-fleld condensation is 
estimated from the decay rate of single scalar fleld (in the vacuum): the energy density of 
the scalar fleld in the condensation is approximated to decrease as 



dt 



'IvPv^ (1-3) 



where 7,^ is the decay rate of (p in the vacuum. However, it has been also pointed out that, 
when the scalar fleld ip is oscillating like Eq. fll.2p . wave functions of flelds which couple 
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to Lp are modified. Consequently, the "decay rate" of tlie scalar field in the condensation 
may be significantly different from the one obtained from the field theory in the vacuum. 
In particular, in some of the cases, instability band may arise in the wave function of the 
final-state particles, which results in catastrophic particle production (so-called parametric 
resonance) [101 [HI [121 [13] . Since the decay of scalar-field condensation is very important, it 
is desirable to have a deep understanding of the decay processes of scalar-field oscillations 
in the framework of well-established quantum field theory. 

In this paper, we consider decay processes of scalar condensations in the framework of 
the quantum field theory. In our analysis, we neglect the effects of cosmic expansion as a first 
step to understand the behavior of the scalar condensation. We postulate that the scalar 
condensation corresponds to the so-called coherent state \ip) in the quantum field theory; 
justification of such a treatment will be also discussed. Then, we will show that the decay 
rate of the coherent state can be systematically calculated by using the unitarity relation 
of the scattering matrix (S'-matrix). By using the fact that the decay rate of the coherent 
state is proportional to the imaginary part of the so-called T-matrix element Q'[((/9|T|<y9)], we 
formulate the calculation of the decay rate of the coherent state. We apply our procedure 
to the case where the scalar-field condensation is coupled to a real scalar field and calculate 
the decay rate for such a case. We will see that our procedure gives the same decay rate 
as that obtained by the discussion of parametric-resonance. We discuss when Eq. ( 11. 3p is 
justified and how the instability band for parametric resonance arises in our framework. We 
also calculate the decay rate of scalar condensation which decays into gauge-boson pair via 
anomaly. 

Organization of this paper is as follows. In Section [21 we first derive basic formulae which 
are used in the calculation of the decay rate of the coherent state. In particular, we define 
the coherent state \^p) in the quantum field theory and present important properties of ly?). 
Then, we explain how the decay rate of the state \ip) is obtained. In Section [31 we consider 
the case where the condensation couples to a real scalar field x via three-point interaction. 
Decay of the coherent state via the interaction induced by the anomaly is discussed in Section 
[H In Section [5l we summarize our results. 

2 Basic Formulae 

In this section, we introduce basic formulae used in our analysis. We discuss how the 
condensation of scalar field decays via some interaction. We assume that the interaction of 
the scalar field is weak enough. 

Total decay rate of any state can be related to the imaginary part of the scattering-matrix 
element due to the S'-matrix unitarity. Let us denote the S'-matrix as 

5 = i + if, (2.1) 

where 1 and T are the unit operator and the so-called T-matrix, respectively. (Here and 
hereafter, the "hat" is used for operators.) Then, from the unitarity of the S'-matrix, the 
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following relation holds: 

f^f ^iif^ -f). (2.2) 

This relation is important for our analysis. 

Wc expect that there exists a quantum state \if) which describes the state with scalar 
condensation. (Details about \(p) will be explained below.) The probability of the state \(p) 
decaying into all the possible final states is related to the imaginary part of the T-matrix 
element as 

Prob(l^) ^ all) = ^ \{f\f\^)\' = 2» [{ip\f\ip)] . (2.3) 
/ 

Now, we consider what the state \ip) is. We first quantize the field operator using the 
free part of the Lagrangian, then treat the interaction terms as perturbations. We denote 
the free part of the Lagrangian of the real scalar field ip as 

jC = ^d.ifd^if - \ml^\ (2.4) 

In our analysis, we use the box normalization of the wave functions with the volume L^. 
Then, the field operator is given by 

- E ^J^, ("xe-^'^^ + cie''^) , (2.5) 



k 



where £'k = y^k^ + m^. The annihilation and creation operators satisfy the following com- 
mutation relations: 

[ak,aL] = 4,k', (2.6) 

while flk and a^/ (as well as and a|.,) commute. 

We postulate that the quantum state describing the scalar-field condensation is the co- 
herent state, which is given by 

|(^) = e-l^^lV2eC^at|Q^^ ^2.7) 

where |0) is the vacuum state satisfying ak|0) = 0. Notice that the state \ip) is properly 
normalized: {ip\ip) — 1. In addition, importantly, the state \ip) is an eigenstate of the 
annihilation operator Sq: 

aoW) = C^\^). (2.8) 

We can also see that 

ip{x) = {ip\^{x)\^) = ip-{x) + ip+{x), (2.9) 
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where 

in r.rl = , ^^^^^ _ ^ 

with 



= -A^e-'^-\ ^ -A;e*--*, (2.10) 



A^ = C^,l ^. (2.111 



One can easily see that, for the coherent state lyj), the expectation value of the field operator 
follows the trajectory of the solution to the classical wave equation. Thus, we expect that 
the coherent state \ip) corresponds to the quantum state where the scalar field is under 
oscillation. 

In calculating T-matrix elements, it is necessary to calculate the expectation values of 
time-ordered products of field operators. By using the Wick's theorem, such products are 
calculated as 

(v^irj]^ (^(xj)|(y9) = ((y9|A^]^ <^(a;j) + (all the possible contractions) 

i i 

= J]^ V9(xi) + (all the possible contractions), (2.12) 

i 

where the symbol T here denotes the time-ordering while is for normal-ordering. In 
addition, in the second equality, we have used Eq. (12. 8p . Even in more complicated cases, 
we obtain 



dip'' 

= Ylfi ifi^i)) + (^11 possible contractions). (2.13) 

i 

Here, we expand fi{ip{xi)) around 0{xi) = (p{xi). Then, {(p\{ip{x) — ip{x))\ip) = 0, and Eq. 
(I2.13P is applicable even when the function fi{(f ) is singular at (f = 0. 

In the following, we will not consider the processes in which (p is produced due to the 
decay of the coherent state. (The inclusion of such processes is straightforward.) In such a 
case, propagator of ip does not show up in the calculation and the field operator (y5(x) can be 
simply replaced by the expectation value (p{x), in which the contraction terms are irrelevant. 

It is also instructive to calculate the expectation values of energy-density operator as 

= L-^ip\Y,E^didM = I^IM, (2.14) 

k ^ 
while the expectation value of the number density is also given by 

= L~^{(f\^dl,dk\ip) = ]-m^\A^\'^. (2.15) 
k ^ 
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For the complex scalar field (which we denote as 0), similar argument applies. We define 
the field operator for the complex scalar field as 

^(^) = E ^J^^ + bW''^) ' (2-16) 

where and 6k (aj, and blj are annihilation (creation) operators. The coherent state for 
the complex field is given by 

10) = e-(l^^l'+l^^l')/2e^^'^o+^^^o|o), (2.17) 

and 

0(x) = (0|0(x)|0) = A^e-'^^' + Ay^^\ (2.18) 

where 

^^ = ^^=, A^ = ^^=. (2.19) 



Energy density of this state is given by 

p<^ = 2m2(|A^|2+|A^|2), (2.20) 

while we can also calculate the expectation values of the number densities of particle and 
anti-particle as 



= L-3(0|^aj,ak|0) = 2m^|Aj2^ (2.21) 

k 

^ L-3(0|^6j,6u|0)=2m^|A^|2. (2.22) 

k 

We can see that the number densities of and its anti-particle are proportional to |^</,P and 
respectively. Thus, when [A^l > [A^l {\A^\ < \A^\), is more (less) abundant than 
in the condensation. It should be also noted that the function (j){x) given in Eq. (12.18^ gives 
an elliptical trajectory on the complex 0-plane. When A^ = or A^ = 0, the trajectory 
becomes a circle and, when \A^\ = \A^\, the trajectory becomes a straight line. 

3 Decay into Scalar Fields 
3.1 Setup 

First, we consider the simplest case where the scalar field ip couples only to the real scalar 
field X via the interaction 

Ant = -^f^fX^, (3.1) 



with /X being the couphng constant. With this interaction, the decay rate of single particle 
in the vacuum is given by 



1 - 



'4m?,' 

A 



(3.2) 



In this section, with the interaction given in Eq. (13. ip . we discuss how the coherent state 
decays. 

As discussed in the previous section, the decay rate of the coherent state can be related to 
the imaginary part of the diagonal element of the T-matrix {ip\T\ip). Importantly, {ip\T\ip) 
is obtained by calculating loop diagrams in the quantum field theory. 

At the one-loop level, in other words, neglecting the fluctuation of ip, {iplTlip) is expressed 

as 



p=l jr{2p) 



(3.3) 



where the summation over JF^^^) is for all the possible Feynman diagrams with 2p external 
if. With p being fixed, one can find 



jr{2p) 



(2p) 



T 



(fx {ip-{x) + ip+{x)} x{x)x{x) 



-|2p 



(3.4) 



where v?- and show up when the field operator ip is contracted with the creation operator 
in \ip) and the annihilation operator in {ip\, respectively. In order to obtain non- vanishing 
results, we need to pick up same number of <y9_ and ip^. 

The right-hand side of Eq. (13.41) contains contribution from various Feynman diagrams 
because and can be ordered in many ways. In Fig. [H we show a typical diagram 
(after imposing the cut). In our notation, we represent the insertion of ip- by the white 
dot o while ipj^ by the black dot •; non- vanishing diagrams have p black and p white dots. 
Internal lines are x-propagator. 

Expectation value of T for the given diagram JF^^p) (with 2p external ip insertion) is given 

by 



UTS:f 



where Sjr is the symmetry factor and 
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2p 



(2vr)4 A=V 



(3.5) 



(3.6) 



In addition, the momentum flowing on the J-th propagator is given by 



j=i 



(3.7) 
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Cut 
ki . 




Figure 1: An example of the Feynman diagram which gives rise to the imaginary part of {ip\T\ip). 
The white dot o indicates ip- insertion while the black dot • is for (^_(_ insertion. In this example, 
p = 8 and, with the cut shown in the figure, N^p = 3. (Notice that other cuts are also possible with 
this diagram.) 



with 



(m^,0) 



and Sj = 1 and —1 if the J-th tp insertion is and respectively. 

The imaginary part of T-^'^*"' is obtained by cutting two propagators (see Fig. [1]), which 
corresponds to the replacements of those two propagators by the 5-functions (with a relevant 
numerical factor): 



53 



T 



jr(2p) 



2p-l 2p 
»^™x i=l j=i+l 



(3.9) 



where 
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Kkl-mk]-Q Uik'j-^ir'. (3.10) 



(For details, see Appendix [Al) 53[7^-^'^*''] is the contribution from the diagram in which cut 



is on i- and j-th propagators. For ^[Tj^ ], we define the energy flow from one side of the 
cut to the other, which we denote -E^xi with the following non-negative integer: 



I=i+l 



(3.11) 



the energy flow is given by 



(3.12) 
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Cut 




Figure 2: Feynman diagram which gives the leading-order contribution to the decay rate in the 
smah amplitude limit. The white dot o indicates insertion while the black dot • is for 
insertion. 

Notice that Q'[7^-^'^*''] contributes only to the decay rate of the process in which of ip in 
the condensation simultaneously annihilate into two x because we neglect the fluctuations 
of the if field. At the perturbative level, such a decay process is kinematically allowed when 
i?vv > 2m V. 



3.2 Small amplitude limit 

In this and the next subsections, we concentrate on the case where the amplitude of (p is 
small. In this case, the leading-order contribution to the E^p = A^<^m0 mode is from the 
diagram with 2iV^ external ip with and being completely separated by the cut. (See 
Fig. O) Concentrating on such a leading-order diagram, the imaginary part of the T-matrix 
element is given by 



Leading 



n [Ck-m'- 
1=1 

Constraints from the 5-functions give kQ 

n [Ck-iQ^r- 



(2vr) 



7=1 



-1 



|A^<^m^. Thus, 



and hence 



53 



Leading 



32n [{Np - 1)!]4 



2ml 



2N^ 



where, for N^m^ > 2m^, the velocity i3m is given by 



\ 



Ami 

A. 



(3.13) 



(3.14) 



(3.15) 



(3.16) 
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(3.17) 



and hence is given by 



Pn^ "^t yuA 



2Ar, 



Leading 



167r [(AT^ - 1)!]4 2ml 



(3.18) 



The above expression is consistent with the result given in the study of the parametric 
resonance [12] . In addition, the decay rate for the A^^ = 1 mode is related to the decay rate 
of single particle, which is given in Eq. (13.21) . as 



Thus, when the amplitude of (p is small, decay of the coherent state can be treated as the 
decay of individual particles in the condensation, which justifies the conventional treatment 
of the decay processes of scalar condensations. 

We also comment here that the decay rate given in Eq. (13.181) is also derived from the 
tree-level calculation of the XX P^ir creation rate in the external oscillating (f field; ^i^^l^^g 
is equal to the production rate of xX P^^t^ P^r unit volume with total energy of E^^ = N^m^. 
In general, at the leading order of A^, the decay rate of the coherent state is also obtained 
by calculating the tree-level production rate of the final-state particles treating the scalar 
condensation as an external field. If we consider higher order contributions, however, such a 
calculation breaks down; treating the scalar condensation as the external field, denominators 
of some propagators vanish in certain types of diagrams. Notice that, in Eq. (13.91) . such a 
difficulty does not exist. (See also the following discussion.) 

3.3 Small velocity limit 

In the previous subsection, we have calculated leading-order contributions to the decay rates 
of each mode in the small amplitude limit. Calculations of the contributions which are 
higher order in the amplitude are straightforward. In this subsection, we discuss when the 
small-amplitude expansion breaks down, taking the A^^ = 1 mode as an example. 

If we calculate 0(|yl^p^) contributions to the decay rate of such mode, which are from 
diagrams with 2p external insertions, one finds that the imaginary part of the T-matrix 
element is inversely proportional to the powers of jSi in the — *■ limit. First, let us 
derive such a behavior with explicit calculation. 

For the = 1 mode, the most important Feynman diagrams in the —>■ 2m^ limit are 
those in which and insertions are next to each other. (See Fig. [31 ) As we will discuss, 
other types of diagrams with fixed p have less singular behavior when /?i ^ 0. We also note 




Leading ^•pI'P^XX- 



(3.19) 
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Figure 3: Feynman diagrams which give the most singular behavior for the N^p = 1 process in the 
small f3i limit. The white dot o indicates ip- insertion while the black dot • is for insertion. 



here that the diagram in Fig. [3] contributes only to the N^p = 1 mode. (The imaginary part 
of the T-matrix vanishes when N^p = 0.) 

Taking into account the diagram shown in Fig. [3l the T-matrix element becomes 



L^T— 
2p 





1 


2 





(3.20) 



and the imaginary part of T^^Zq is given by 



$5 



"32^""^ 



(4p-7)!! 



fiA, 



2ml(3j 



2{p-l) 



+ 0{(3r'n. (3.21) 



p\{p — 1)! 

where {2p - 1)!! = W=i{21 - 1) for p>l, and (-3)!! = -1. (For details, see Appendix (B) 
Hereafter, we neglect 0{l3l~^^) contribution in Eq. fl3.2ip . and the decay rate of the coherent 
state becomes 

2(p-l) 



IGtt 



(4p-7)!! 



E 



lp!(p- 1)! 



2mlP! 



(3.22) 



As discussed in Appendix IHl the inverse powers of f3i stems from the derivative of the 
function B{rn?^] ^i, ^j) given in Eq. ( 1A.7I) with respect to or ^j. The order of the derivatives 
is equal to the number of the propagators whose denominators vanish in the on-shell limit 
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1 — 2p). Number of such propagators is maximized for the diagram given in 



Fig. [31 Thus, we safely neglect other types of diagrams in studying the case of (3i 0. 



As we mentioned, F^^^q^'' becomes singular when (3i 



0. In other words, for the = 1 



mode, the small-amplitude expansion breaks down when |y4<^| is comparable to fi~^m'^/3f. 
These behaviors are related to the fact that instability bands appear in the solution to 
the classical wave equation of the scalar field (i.e., x i^i our argument) which couples to a 
oscillating scalar field (i.e., ip). With the interaction given in Eq. (13. ip . the wave equation of 
X (with the 3-momentum k) in the oscillating background is given by the Mathieu equation: 



^1 cos nipt 



0. 



(3.23) 
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Parametrizing the momentum of x 



= -ml (Pl + e , (3.24) 



1 

the lowest instabihty band in the small-amphtude hmit is given by 

-e<e<e, (3.25) 

where 

e = (3.26) 

From the study of the parametric resonance, the momentum of x produced by the decay 
of the scalar condensation is in the range given by Eq. (13.251) . For the consistency of the 
calculation, the mass of the initial-state particle (p should be large enough so that is 
positive even for e ~ —6; otherwise, 0{6^) contributions may be also important. This 
argument gives the limitation of the small-amplitude approximation; Pi ^ 9^^"^ is required, 
which results in \A^p\ 3> fi~^mlPf. 

The above argument is supported by the fact that F^^^q^'' given in Eq. fl3.22p is equal 

to the decay rate obtained by the parametric-resonance analysis. Indeed, F^^^q^'' is also 
expressed as 



" 2^^y|k|_<|k|<|ic|,(27r)3^ 



1287r 



^< £ de^{P! + e) (e2_e2), (3.27) 



where, in the first equality, |k|-|- = \m^^ jdl ± 9, and A*^^*" = \^/9'^ — is the "growth- 
rate factor" for the A^^ = 1 mode obtained in the study of the Mathieu equation [T2j. The 



equivalence of Eq. (13.221) and Eq. (13.271) can be seen by expanding y /J^ + e in the integrand 
of Eq. (I3.27P around e = (assuming jSi > 9^/"^). Eq. (I3.27P is nothing but the decay rate of 
the scalar condensation in the small-amplitude limit derived from the parametric-resonance 
analysis. 

Eq. (I3.27P (and Eq. fl3.18p ) also shows the fact that, at least at the small-amplitude limit, 
the results from the parametric-resonance analysis is obtained in our procedure where the 
quantum state describing the scalar condensation is postulated to be the coherent state. The 
equivalence of two approaches is also expected from the fact that the basic equations gov- 
erning the behavior of the parametric resonance is derived in our framework. In particular, 
we can calculate the density matrix of the final-state particle x the quantum field theory. 
We can see that the Mathieu equation shows up in the calculation and that the resultant 
density matrix is the same as the one obtained in the study of parametric resonance. These 
subjects will be discussed in the next subsection. 
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3.4 Calculation of the density matrix 



In the classical treatment of the parametric resonant system, it is well known that the 
Mathieu equation appears as the equation of motion for the x field. Thus, the equation is 
also expected to be obtained in the quantum-field-theory treatment. In this subsection, we 
show the derivations of the Mathieu equation and the density matrix of x explicitly. 

Since we are interested in the case where the scalar field initially forms the scalar 
condensation oscillating around the minimum of its potential, we describe the initial state 
(which is taken at t = in this subsection) as 

V) = \^) ® |0)x, (3.28) 

where the first and second kets represent the states for (f and Xy respectively. In addition, 
|0);^ is the vacuum of the x field. (In the following, we omit the subscript x-) 

For our argument, it is convenient to use the density matrix of the total system in the 
Schroedinger picture. The density matrix at the time T is simply given by 

Aot(r) = e-^^^K)(^|e^^^, (3.29) 

with H being the Hamiltonian of the total system. We consider the properties of ptot in the 
coordinate basis: 



\q(^X) = \q)(E)\X), (3.30) 

where \q) and |X) are eigenstates of the field operators and x, respectively: 

(^(t,x)|g) = g(x)|g), x(t, x)|X) = X(x)|X). (3.31) 

Then, the density matrix in the coordinate basis, ptot[Q, X; q', X'] = {q ® X\ptot(T)\q' ® X'), 
is given by 

pUq,X;q',X'] = JVq.J Vq[ j VX, J VX'M^) {X,\0) {^\q';) {0\X[) 

K[q,X;q,,X,]K*[q',X';q'„X'^. (3.32) 

The kernel is represented in the path integral form as 

K[q,X;q,,X,]= Vx Py^e*^-*, (3.33) 

./x(0,x)=X,(x) i</.(0,x)=g,(x) 



where S'tot is the total action. For the explicit form of the kernel, see [TH [15] . 

In order to study the behavior of the x field, we derive the reduced density matrix of x 
by tracing out q and q' variables: 

p,ed[X; X'] = j Vqptoth X; g, X']. (3.34) 
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For this purpose, we define 



(3.35) 
(3.36) 



Then, with the use of the properties of the coherent state, we obtain the following relations: 
{q\(p) = exp -l-m^{qo - L^^^A^)^ Yi exp -^^kl^kP , (3.37) 

(X|0) = nexp[-^^k|Xkpj, (3.38) 

where, in this subsection, we consider the case where is real. (A^ can be taken to be 
real with a relevant shift of the time variable.) In addition, = + m^. In Eqs. (13.371) 
and (13.381) . we omit unimportant numerical constants. Using Eq. (I3.37p . the reduced density 
matrix becomes 



Pred[X;X'] 

where 



vx,vxl{x,\o){o\x:) 



x(r,x)=x 
x(o,x)=x. 



x'(T,x)=X' 
X'(0,x)=x; 



Vx'e 



'JS[x]-iSlx']+C[x,x'] 



Six] 



dt / d 



^df^xd^X - ^^Ix^ - ^fiA^X^ cos m^t 



,(3.39) 



(3.40) 



while C gives the collision terms: 

.2 .T 



c[x,x'] 



" 4 Jo 

X^{x) 



dt / dt' / rf^x / d 



d^k 



(27r)32Ek 



cos k (x — x') 



X'\x) 



(3.41) 



It can be seen that the periodic perturbation term, x^ cos m<^t, appears in the reduced density 
matrix. Collision terms, which are proportional to x'^i ^ind ^ilso appear after the 

integration. Since no approximation was made to derive the density matrix, the above 
formula can be used at any amplitude of the field, A^. Furthermore, we can derive the 
kinetic equation by calculating the correlation function of the x fi^ld on the reduced density 
matrix, which allows us to describe the non-linear dynamics of the x system, as pointed out 
in [IE]. 

When the amplitude A^p is small enough, we can treat the collision terms as perturbations. 
Then, the reduced density matrix is, at the leading order calculation, written by using the 
wave functional of the x field: 



Pred[X;X'] = vl/[T,X]xvl/*[r;X'] 



(3.42) 
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where 



x(T,x)=X 
x(0,x)=X, 



Vxe 



(3.43) 



The wave functional is written as a product of an infinite set of wave functions of harmonic 
oscillators as 



-uj]^\X, 



i,k| 



Xk(0)=Xi,k 



where 



dt 



1,^ ,2 



(3.44) 



(3.45) 



with the "dot" being the derivative with respect to time. \E'[T, X] satisfies the boundary 
condition \1/[0,X] = Hk ^"'^'''"^'''^^^ (up to normalization), and its evolution is governed by 
the wave equation derived from the action given in Eq. (I3.45p . Thus, "^[T, X] is nothing but 
the wave functional obtained in 



k \/Uk{T) 



exp 



where U]f_{t) is the solution to the Mathieu equation: 



(3.46) 



+ (ul + fiA^ cos m^tj Mk = 0, 



(3.47) 



with the conditions Mk(0) = yTr/uj]^ and Uk(0) = iy/uj^- This fact supports that the scalar 
condensation is well described by the coherent state in the quantum field theory. 

Before closing this section, we emphasize that the quantity Q=[(v9|T|(y9)] is relatively easily 
calculated with wide variety of interactions and final-states. Thus, for some applications, 
our procedure is more powerful than the approach using the Mathieu equation. 



4 Decay via Anomaly 

Next, let us consider the case where a complex scalar field may decay via chiral and 
conformal anomalies. As the fundamental theory, we expect that there exists chiral fermions 
which have gauge quantum numbers and that the complex scalar field couples to the chiral 
fermions through a Yukawa interaction. To make our discussion definite, we consider SU (Nc) 
gauge interaction; chiral fermions Ql and are in fundamental and anti-fundamental 
representations of SU{Nc), respectively, while is singlet. 

The complex scalar field couples to chiral fermions Ql and via the Yukawa interaction 

vCvukawa = -|/(0(5lQr + h.c). (4.1) 
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It should be noted that, in this model, there exists anomalous U{1) symmetry, which we call 
U{1)a', charges of 0, Ql and are 1, — |, and — |, respectively. 

When the amplitude of (p is large, fermions Ql and acquire Dirac mass. Thus, when 
"^(j) ^ 2/101) effective mass of the fermions are much larger than the mass of (p. In this case, 
the decay process — > QlQr is expected to be kinematically forbidden. 

When m,^ <^ y\4'\, it is rather convenient to consider the low-energy effective field theory 
by integrating out the fermions. The relevant (light) fields in the low-energy effective field 
theory are (f) and gauge fields as far as the effective mass of the fermions are much larger than 
m^. In the following, we concentrate on such a case; thus, we assume that the inequality 

-C y\(j){x)\ always holds at any point of the trajectory of (p. 

We first consider the effects of the operator induced by the chiral anomaly: 



£efr = -^Ai (ln0 - ln0tj p^'^F^,, (4.2) 
where F^iy is the field-strength tensor and 

^ leM^p-iT^^. (4.3) 

In Eq. (14. 2 p and hereafter, summation over the adjoint gauge index is implicit. In addition, 

A, ^ J^r«. (4.4) 

with g being the gauge coupling constant of SU{Nc), and Tr = |. 

Now, let us discuss the decay of the coherent state given in Eq. fl2.17p . At the leading 
order in Aj, which is of 0(Aj) in the calculation of (0|T|</)), we obtain 

(0|2r|0) = --XlJ d^xd'x'{0\TF,,{x)F {x)F,,,,{x')F {x')\0) 

\n (j){x) — \n (j)'^ (x) ln0(x') — ln0^(x') . (4-5) 
In the following, we consider the case that > A^. Then, we expand ln0(x) as 

ln(f){x) = ln(A^e-''3^") + ^ ^ e^™^^", (4.6) 

n=l ^ \^<t>/ 

where = {m^,0). At the lowest order in Aj, the decay rate of the coherent state can be 
obtained by calculating the two-point functions of several types of operators with relevant 
momentum injection. For the local operator 0{x), let us define 

Xo(Q) = -iJ d^Xid^X2{0\Td{xi)d{x2)\0)e''^^"'\ (4.7) 
Then, by using the fact that F^^F^" is expressed as a total derivative: 



£^^"^,(9^^,) + (gauge field)=^ , (4.8) 
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Eq. (14.51) becomes 



-2Ai%.^^(0). 



oo 
n=l 



A, 



2n 



I 



(4.9) 



It is notable that Xqmj^^(O) has no imaginary part because there is no momentum injection 
into the internal gauge-boson lines from the Q^-ft'^-vertex. Thus, the coherent state does 
not decay if = or = 0. (This statement holds even after taking into account the 
higher order terms in Aj.) This fact can be understood by the conservation of the U{1)a 
charge. With a fixed value of the total energy of the system, U{1)a charge is maximized when 
A^ = Q OT A^ = 0. Thus, if U{1)a charge is conserved, the decay of in the condensation 
into the gauge bosons is forbidden. Of course, the interaction given in Eq. (14.21) breaks U{1)a 
symmetry because £efr is not invariant under the U{1)a transformation. This is due to the 
fact that £eff is induced by the chiral anomaly. However, we can add new fermions, which 
we call Q'l^ and Q'^, to have conserved U{1) symmetry. Indeed, with and Q'^, which are 
in fundamental and anti-fundamental representation of SU{Nc), respectively, we can define 
non-anomalous U{1)a symmetry by assigning charge +i to both of Q'^^ and Q'^. (Notice that 
Q']^ and Q'^ do not have to couple to 0.) In this case, conservation of the U{1)a charge is 
obvious and the decay of into the gauge bosons is completely forbidden. Thus, the U{1)a 
charge stored in the scalar condensation cannot be released by the interaction given in Eq. 
(14. 2p . This fact may have some relevance in the study of the decay of scalar condensations in 
various cosmological scenarios, in particular, in the Affleck-Dine scenario [9j . In the absence 
of Q'l and Q'^, instanton effects may generate new interactions which explicitly breaks U{1)a 
symmetry. In such a case, decay of the coherent state occurs via such new interactions. 

Using the relation 



53 



'-^F,,^Ft^'^{Q) 



An 



the decay rate is given by 

r 



FF 



^ n=l 



n 



A, 



A, 



2n 



(4.10) 



(4.11) 



Since the decay rate vanishes if A^ = or A^ = 0, the decay of the coherent state in this 
case should be understood as an annihilation between (p and (p in the condensation; same 
number of and annihilate into the gauge boson pair. 

In the study of the decay of coherent state, the energy-loss rate is also important. Using 
the fact that the imaginary part of 2p^^p^r,{2nQ^) is from the decay process into two gauge 
bosons with the total energy of 2nm^, the energy- loss rate can be calculated. So far, we 
have considered the case where A^ > A^. However, the decay rate for the case of A^ < A^ 
is derived by interchanging A^ ^ Ai in the result. Thus, we obtain 



dt 
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FF 



71 



(iV2 - l)Xlm 



oo 
ra=l 



(4.12) 
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where we have used the fact that the number densities of and are proportional to I^^P 
and respectively. (See Eqs. (JMH) and (^2^.) One may simplify the above energy-loss 

rate by using 



oo 

E 

n=l 



r(l + 4r + r^) 
(1 -r)4 



(4.13) 



Here, we emphasize that the results given in Eqs. (14. lip and (I4.12p can be used for any 
value of the amplitude (as far as the effective mass of Ql and are much larger than m^). 
When ^ or <^ n^, which corresponds to the case where the classical motion of 
the scalar condensation is almost circular, the energy-loss rate is well approximated by the 
leading term in Eq. (14.121) . On the contrary, in the limit of — n^, higher order terms 
become important and the energy-loss rate is enhanced. In this case, however, one should 
note that, at some point of the classical trajectory, |0| approaches to the origin. Then, the 
effective mass of the fermions Ql and may become so small that the effective field theory, 
which is obtained by integrating out these fermions, may break down. We also note that, 
with the ratio n^/n^ being fixed, the decay and energy- loss rates are independent of the 
amplitude of the scalar condensation. 

Before closing this section, we also present the result for the case where the scalar field 
(p couples to the gauge field as 



£eff = AR (ln0 + ln0t) F^^F^,. 



(4.14) 



This type of interaction is also generated by integrating out particles which acquire masses 
from the condensation of (f) (like Ql and Q%)- At the leading order in Ar, the energy- loss 
rate is given by 



dt 



FF 



'An! 

71 



oo 

E 

n=l 



n 



max(n<^, n^) 



(4.15) 



5 Summary 

In this paper, we have discussed the decay processes of the scalar condensation. We postu- 
lated that the quantum state corresponding to the scalar oscillation is the so-called coherent 
state in the quantum field theory. Then, by using the 5'-matrix unitarity, we have developed 
the method to calculate the decay rate of the coherent state. We believe that our procedure 
can be applied to a large class of models which may contain various types of interactions. 

Then, in order to demonstrate how the decay rate is calculated, we considered two exam- 
ples. First, we studied the case where the scalar field ip couples to another scalar field x via 
three-point interaction. Using the small-amplitude approximation, we have calculated the 
decay rate for the process where {N^ = 1, 2, 3, ■ ■ ■) of (p in the condensation simultane- 
ously annihilate into a pair of x- For the case of = 1, we have seen that the result is the 
same as that in the conventional approach where the decay rate of the scalar condensation 
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is estimated by the product of the decay rate of single if in the vacuum and the number 
density of ip. We have also pointed out that the small- amplitude approximation breaks down 
when the amplitude becomes close to fi~^m'^Pf, where is the couphng constant and j3i is 
the velocity of x iii the iV^ = 1 mode. Such a behavior is also expected from the discussion 
based on the parametric-resonance. Indeed, our procedure reproduced the decay rate of the 
scalar condensation calculated from the parametric-resonance analysis. 

The second example was the case where the complex scalar field decays into gauge bosons 
via the interaction induced by the chiral anomaly. We have considered the case where the 
scalar potential has U{1)a symmetry to rotate — >■ e*"0 at the classical level, which is 
broken by the effect of the chiral anomaly. In this case, we could calculate the decay rate 
without using the small-amplitude approximation. We have seen that the decay process is 
forbidden unless both the particle and its anti-particle exist in the condensation, and 
that the "decay" of the coherent state is due to the annihilation between them. Thus, U{1)a 
charge stored in the condensation cannot be released by the effective interaction induced by 
the chiral anomaly. 

In our analysis, the effects of the cosmic expansion were completely neglected. How- 
ever, we believe that our results are applicable to the cosmological discussion as far as the 
expansion rate of the universe is smaller than the mass of the scalar condensation. 

Acknowledgement: This work was supported in part by the Grant-in- Aid for Scientific Re- 
search from the Ministry of Education, Science, Sports, and Culture of Japan, No. 19540255 



A Calculation of the Imaginary Part 

Although the technique which will be explained here is well-known, in this appendix, we 
show how Eq. (13. 9p is derived for the sake of some of the readers. For this purpose, we 
calculate the imaginary part of the following quantity: 



(T.M.). 




(A.l) 



Here, 



I 



k + J2 ^jQ' 



(A.2) 



j=i 



where Si = ±1 (with -|- ■ ■ ■ -|- = 0), and = (m<^, 0). 

In order to calculate the imaginary part of /jr, it is convenient to rewrite /jc- as 




(A.3) 
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where the hmit ^ indicates that C,i —* to^ (J = 1 — 2p); before taking the hmit, ^/ are 



all set to be different. Integrand of Eq. flA.3l) has poles at /cq = ~"^ipZ]j=i^j i yk/ + 
(/ = 1 — 2p) and, after fco-integration, Ijr becomes 



2p 

-TT lim ^ / 



X i=l 
2p 



2k 



/^i 



-vr 



(A.4) 



where, in the second equality, the fco-iiitegration is performed in the region where kio > 0. 
Using the relation {x + iO~^)~^ = P{x~^) — iTT6{x) (where "P" is for the principal value), we 
obtain 



2p-l 2p 



-X i=l j=i+l ' 



(A.5) 



Substituting the above expression into Eq. fl3.5p . we obtain Eq. (13.91) . 



With the quantity 



ki and fcj are related as ki = kj + N^Q^ (or 
/cj = kj — N^Q^). Thus, if N^p = 0, ki = kj and the imaginary part vanishes. Notice also 
that the constraints from the (5-functions can be solved; by shifting the integration variable 
kj ki and kj can be taken to be k and k — N^Q^, respectively. Then, the constraints from the 
(5-functions become kQ^ = ^Ny^m'^. Consequently, the product Hj-^i j^kj — becomes 
fc- independent. The remaining part is proportional to the two-body phase space for the 
process where the parent particle with mass A^^m^ decays into two daughter particles with 

^1/2 , ^1/2 

masses 4i and 4 : 



327r3 



5(ArJmJ;e„e,), 



(A.6) 



where, for ^/Q^ > ^i + ^ 



1/2 >.l/2 



3 ' 



(A.7) 



while B{Q'^;^i,^j) = for \^Q^ < + ^j^^. Notice that the function B is related to (3^^ 
given in Eq. fl3.16p as 



PN=B{N'm';m'm'). 



(A.8) 
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B Derivation of Eq. ( 13:211 ) 

In this Appendix, we calculate the imaginary part of the following integral 



(P -ml + zO+) ' [(k - Q^f -ml + iQ' 



(B.l) 



to derive Eq. fl3.2ip from Eq. fl3.20p . Using the procedure given in Appendix |X1 we express 
^FigEjas 



—I 



€'''^"'x- \^"J 1=1 J=l 

The imaginary part of this quantity is obtained with 



(B.2) 



1=1 



-1 



i=l Ij^i 



(B.3) 



and with the similar replacement of the second product. Then, the imaginary part of -^Fig[3] 
becomes 



FigE] 



p p 



We can use the relation: 



lim ■ ■ ■ lim E/(^i) IT' 



dP- 



3? ]_ '^'P 



(p - 1)! rfxP- 



(B.5) 



to obtain 



Figl3] 



167r [(p-1)!] 



(B.6) 



Taking the 0(/3[ ^'^^ term from the above expression, which is the most singular one 
when A ^ 0, Eq. (IX^ is derived. 
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